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Abstract 

An Important consideration In the global local finite-element method 
(GLFEM) Is the availability of global functions for the given problem. 

The role and mathematical requirements of these global functions In a GLFEM 
analysis of localized stress states In prismatic structures are discussed. 

A method is described for determining these global functions. Underlying 
this method are theorems due to Toupln and Knowles on strain energy decay 
rates, which are related to a quantitative expression of Sal nt- Vena nt' s 
principle. It Is mentioned that a mathematically complete set of global 
functions can be generated, so that any arbitrary Interface condition 
between the finite element and global subregions can be represented. 
Convergence to the true behavior can be achieved with Increasing global 
functions and finite-element degrees of freedom. Specific attention Is 
devoted to mathematically two-dimensional and three-dimensional prismatic 
structures. Comments are offered on the GLFEM analysis of NASA flat panel 
with a discontinuous stiffener. Methods for determining global functions 
for other effects are also Indicated, such as steady-state dynamics and 
bodies under Initial stress. 
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Introductl on 


The finite-element method (FEM) has revolutionized structural and 
stress analyses In the last quarter of this century. Its Impact has been 
widespread, even extending beyond the preserve of structural engineers to 
other fields. Although FEM Is acknowledged as an extremely powerful 
modeling technique, the analysis community with Its collective experience 
will admit that It Is not the quintessential technique. There are problems 
not well suited to FEM that result In clunsy, Ineffective and costly 
mathematical models. Examples can be cited from problems involving stress 
singularities and Infinite domains. To obviate the difficulties, 
modifications to FEM have been explored. One alternate approach which 
bodes considerable promise Is the so-called Global-Local Finite-Element 
Method (GLFEM) . 

GLFEM utilizes both conventional finite elements and classical Rltz 
functions In the modeling process. Their respective roles are readily 
apparent; finite elements work well In regions where complicated geometry 
and Inhomogeneous material characterizations prevail, and Rltz functions, 
hereinafter referred to as global functions In GLFEM, enable the behavior 
other regions to be represented accurately and efficiently. At this stage 
of development, GLFEM can be assessed to be In Its maturing phase. It Is 
of good lineage, has already exhibited an enhanced capability above FEM In 
certain problems, and promises effectiveness In other classes of problems 
upon Its full development. 

Herein, GLFEM as applied to the analysis of localized stresses In 
prismatic structures Is discussed. First, the essence of GLFEM and 
various GLFEM modeling layouts are sunmarlzed. A brief review of some 
problems that have been successfully analyzed by GLFEM Is given. Then, 
the main theme relating to GLFEM analysis of localized stress states Is 
addressed. Prismatic structures that can be described mathematically 
by two spatial variables are discussed first. Attention Is devoted to 
the global functions, their development and their roles In the present 
setting. Then, three-dimensional structures are considered, with 
reference to the NASA example problem, where an outline of a method of 
attack Is given. Last, comments on the analysis of localized stresses 
Involving steady-state dynamic effects as well as other conditions are 
given. 


Basic Concepts of GLFEM and the Various Mesh Configurations 

Hamilton's principle, or alternatively the theorem of mlnlmun 
potential energy when no Inertial effects are present, may be considered 
as the basis for generating GLFEM equations. The theory and variational 
derivation of these equations may be found In Ref. [1, pp. 451-474] . Also 
Included therein Is a survey of GLFEM contributions to the literature up 
to 1982. 

As noted earlier, the technique utilizes finite-element modeling with 
classical Rltz approximations simultaneously. It enjoys the advantages of 
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more versatile modeling capabilities with substantially fewer degrees of 
freedom. Various global /local modeling configurations are Illustrated In 
Fig. 1. Figs, la and If represent* respectively* the classical Rltz and 
finite-element configurations. The others are possible GLFEM mesh layouts. 
In a given problem, the modeling may take the form of any one of these con- 
figurations or a combination of two or more of them for various subregions. 
An Important key Is the enforcement of kinematic Inter- regional continuity 
between various global and local subregions by means of constraint equa- 
tions. In problems on localized stress states* only the Fig. lc configura- 
tion will be used* where finite elements exclusively are used In one sub- 
region and global functions In the other. Moreover* the global subregion 
may be Infinite In extent. 

The governing matrix equations In a GLFEM analysis have the form: 
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where {6} denotes the finite-element degrees of freedom and (SI contains 
the array of generalized coordinates associated with the global functions. 
In Eq. (1), CKgg] * [Mgg3* CK tl 3 » and CM^ 3 refer to the global and local 
stiffness and mass matrices of the system. The matrices [K^g] = CKj^g]' 
and CMjtg 3 = CM^g]^ represent global -local coupling from Imposing kinematic 
continuity at Interfaced ) between subreg1on(s ). Details on the formation 
of these matrices may be found In Ref. Cl], 

It Is mentioned that GLFEM variants are possible* which do not lead 
to the same set of governing equations as Eq. (1). These variants contain 
the spirit of GLFEM and employ the modeling configurations shown In Fig. 1; 
however, the method of enforcing 1 nter- regional contl nulty may differ. An 
application concerned with elastic wave scattering will Illustrate one such 
varl ant. 

Another key point In GLFEM Is the availability of an appropriate set 
of global functions for a given problem or a class of problems. The 
accuracy and effectiveness of the method are dependent upon the quality 
of the global functions. The choice of these global functions for the 
analysis of localized stresses In prismatic structures and their method 
of derivation will be discussed In what follows. It will become apparent 
why these global functions* together with the finite-element model of 
the subregion that contains the localized stresses* will lead to a 
superior model. 


Seme Examples of Global Functions for GLFEM 

Two areas Ideally suited to GLFEM are fracture mechanics and Infinite 
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and/or semi-infinite domain problems. Much has been published on various 
aspects of fracture mechanics problems. Many numerical methods have been 
used, many falling within a GLFEM classification or Its variants. The 
global subregion model usually takes the form of special crack tip 
elements* where the singular stress field is incorporated Into the 
stiffness matrix. These elements are well known. No further elaboration 
on this subject will be given here. Elastostatlc analyses of half-space 
problems, wherein the far field behavior Is represented by global functions 
(for example, Boussl nesq or Cerruti solutions), have also met with 
considerable success. A number of references on both of these subjects 
may be found In Ref. [1], 

Herein, two recent GLFEM applications are mentioned to emphasize the 
roles of the global functions and their mathematical suitability. They are 
concerned with (1) steady-state elastic wave scattering by axlsymmetrlc 
objects embedded In an Infinite isotropic medium and (2) steady-state 
soil -structure Interaction Involving an axlsymmetrlc structure occupying 
some locale In a semi-infinite medlun. The feature of note Is that these 
global functions constitute a complete set of eigenfunctions and have the 
capability of mathematically representing an arbitrary scattered field to 
any given accuracy. Hence, the true behavior In the far field can be 
achieved. 

In Fig. 2 Is shown an elastic, axlsymmetrlc Inclusion embedded 1 n an 
elastic. Isotropic medlun. Because finite elements are used for the object 
It may have Inhomogeneous, orthotropic properties. The finite-element 
subregion Includes this object and a portion of the surrounding medium. 

For convenience In the analysis, the Interface Is taken to be spherical. 
Outside of the finite-element subregion Is the outer field, where a 
complete set of outgoing spherical harmonics Is used to model the scattered 
field. Each component satisfies the equations of motion and the Scmmerf el d 
radiation conditions. The global functions have specific stress and 
displacement distributions at the Interface, and their undetermined 
strengths are the global function coefficients or the terms In {SI. A 
given Incident wave 111 uni nates this object. The scattered field Is 
determined by solving the finite-element equations and requiring that the 
sun of Incident and scattered wave fields based on the global functions 
have both traction and displacement continuity with the finite-element 
data at the Interface. Details of this analysis may be found In Ref. [2]. 
Here, attention Is called to the mathematical flexibility of the global 
functions for accommodating Interface continuity to any precision with a 
sufficient number of terms. 

The dynamic soil -structure interation problem under steady-state 
conditions Is shown in Fig. 3. The approach used here Is similar to that 
for elastic wave scattering by an object embedded In the entire space. 

In fact, the same set of spherical harmonics for the entire space may be 
applied to this half-space problem. However, traction-free surface 
conditions are not satisfied by the spherical harmonics. Thus, In 
addition to traction and displacement continuity at the hemispherical 
Interface, It Is necessary to enforce the traction-free surface In the 
global subregion. In Refs. [3,4], details concerning an Integral 
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constraint condition to meet this traction-free surface condition are 
given. Again, It Is noted that because a complete set of global 
functions Is used (that Is, a set capable of modeling any arbitrary traction 
and displacement conditions between various subregions In a GLFEM layout), 
the analysis procedure enjoys the opportunity of converging onto the 
true behavior with Increasing FEM and global degrees of freedom. 


Mathematically Two-Dimensional Structures 

The choice of global functions for mathematically two-dimensional 
structures will now be discussed. As Illustrations of this class of 
problems and their GLFEM layouts, refer to Fig. 4, where examples of a 
laminated composite plate and cylinders are given. The double lap joint 
may be considered as a plane strain problem herein. The scarf joint 
joining two cylinders may be taken as an axl symmetric structure under 
axl symmetric or asymmetric loads. The purpose Is to study the stresses 
In these joints. 

Uniform stress states exist at points well away from these localized 
stress regions. If FEM were used. It Is obvious that an awkward model 
would result. In GLFEM, two-dimensional finite elements (planar or 
axl symmetric toroidal elements) are used for the subregion containing the 
localized stresses. If the localized stress state contains a singularity, 
a global subregion within the finite-element subregion may be added. The 
Interface location Is dependent on the global functions' mathematical 
capability for capturing the transitional stress and displacement fields 
accurately. For global functions capable of representing the true 
behavior, the finite-element subregion can be quite small with the 
Interface(s) near to the localized stress area. An Independent set of 
global functions must be adopted at each Interface. For the lap joint 
In Fig. 4, two or three distinct systems of global functions may be 
needed depending on the thickness and material properties of the plate 
components. Each set of global functions Is associated with Its own 
set of generalized coordinates or global coefficients. For the 
cylindrical scarf joint, two Independent sets are needed. 

The global functions In these cases are based on theorems relating to 
a quantitative expression of St. Venant's principle. Toupln [5] and 
Knowles [6] presented upper bound estimates of strain energy decay rates 
In terms of distance from a self-equilibrated stress state. Their 
results can be stated In the form of a strain energy Inequality: 

V (x) « V(0) e‘ 2Yx (2) 


where Y Is the Inverse of the characteristic decay length, V (0 ) Is the 
total strain energy and V(x) Is that portion of V (0 ) In the body beyond 
x. Since the strain energy Is quadratic, the mechanical variables such 
as stress, strain and displacement are of the forms: 
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Uj(x) < K^e“ YX e^.fx) * K 2 e~ YX a^(x) * K 3 e_YX (3) 
where (1=1»2»3) are constants. 

Based on these theorems* a boundary-value problem can be formulated 
for a prismatic structure. Using the solution form In the prismatic 
direction as e"Y x , the analysis leads to an eigenvalue problem. The 
eigenvalues y's are the characteristic Inverse decay lengths and the 
eigenfunctions are distributions of sel f-equll Ibrated stress states. 

These elgendata comprise a complete set from which any arbitrary self- 
equilibrated stress state may be represented. These elgendata may be 
used as global functions for describing the far-fleld behavior In a 
prismatic structure. Horgan and his colleagues have solved a number 
of problems on homogeneous and sandwich plates under plane strain using 
the Airy stress function as the primary dependent variable (see* for 
examples* Refs. [7*8]). 

For a laminated composite structure* It Is more convenient to 
determine the elgendata numerically. Dong and Goetschel [9] developed 
a one-dimensional finite-element analysis for extracting elgendata for 
a laminated composite plate with an arbitrary number of bonded, elastic 
laminates. Finite-element discretization occurs in the thickness 
direction* see Fig. 5. Applying the theorem of minimum potential energy* 
a system of second-order ordinary differential equations Is obtained. 

By Invoking exponential decay e' x , the following second-order algebraic 
eigenvalue problem results: 


[K-j ]{Q) - yCK 2 ]{Q} + y 2 [K 3 ]{Q} = 0 (4) 

where {Q> Is an ordered set of the plate’s nodal displacements. This 
equation Is reducible to first order with a non-symmetrlc matrix. If a 
large number of degrees of freedom are Involved, a Block-Stodol a Iteration 
technique [10] can be used to extract the elgendata efficiently. The 
solution consists of a complete set of eigenvalues and corresponding 
eigenvectors, which are the self-equilibrated displacement states for the 
given composite plate. Stresses can be computed from these displacements. 

Laminated cylinders may also be solved using the same finite-element 
scheme* see Ref. [11]. The mechanical variables have circumferential 
dependence* which may be expressed analytically by Fourier series. As 
a circumferential mode number m occurs In this case* the counterpart to 
Eq. (4) for each circumferential mode has the form: 

[K <m)]{Q} - Y[K 2 (m)]{Q} + Y 2 D< 3 (m)]{Q} = 0 (5) 

The solution to Eq. (4) or (5) provides the global function data 
base for the numerical evaluation of the global stiffness matrix and the 
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global-local coupling as a prelude to a mathematically two-dimensional 
GLFEM analysis of the localized stress zones. Sane preliminary results 
of this type have been obtained, which are contained In Refs. [12,13]. 

These limited scope studies Indicate an overall feasibility for this 
approach. 

The accuracy of the global functions depends on the fineness of the 
one-dimensional finite-element model adopted for the elgenproblem. Since 
only one-dimensional finite elements are used, a large model does not Incur 
an Inordinate computational effort because of a very small bandwidth. 

The number of global functions required In a GLFEM analysis depends 
on both the nature of the localized stress and the location of the 
Interface. Having an Interface near the localized stress zone will 
require a larger number of global functions, but with a decrease In the 
finite-element coordinates. Conversely, an Interface far removed from 
the localized zone needs fewer global functions, but Is counteracted by 
a greater nunber of finite-element degrees of freedom. 


Three-Dimensional Structures and the NASA Problem 

A schematic of a three-dimensional prismatic structure and the NASA 
problem of a flat stiffened composite panel with a discontinuous stiffener 
are shown In Fig. 6. In this class of problems, three-dimensional finite 
elements must be employed In the localized stress region. Global functions 
must be used at the Interface. They can be obtained from a two-dimensional 
finite-element analysis of the Inverse characteristic decay lengths. 

The analysis to determine the global functions follows the same 
methodology as that for mathematically two-dimensional structures. The 
prismatic cross-section Is modeled by two-dimensional finite elements. 

With the dependence In the prismatic direction taken as e"Y x , a n 
elgenproblem emerges for the extraction of eigendata that form the 
global function data base for the given cross section. The other 
aspects are the same as that described In the previous section. It Is 
obvious that, In this case, the computational effort Is greater. 

Some comments can be given on a GLFEM analysis of the NASA flat panel. 
The set of two-dimensional global functions constitutes a complete system 
of eigenfunctions, with the non-zero eigenvalues associated with Inverse 
characteristic decay lengths of self-equilibrated stress states. There 
are two zero eigenvalues for two stress distributions exhibiting no decay. 
They are the uniform axial deformation and pure bending states. These two 
global functions are needed In a GLFEM analysis of the NASA flat panel, 
since the discontinuous stiffener may produce bending In addition to Its 
uniform end shortening. A set of global functions with all of these 
members present should permit a three-dimensional finite-element model 
to be concentrated on the details of the discontinuous stiffener region. 
Parametric studies wherein the hole and the gap length In the 
discontinuous stiffener are varied may be conducted. Each configuration 
will require a change of the three-dimensional finite-element mesh, but 
the same set of global functions may be used in all cases. 


Applications to Steady-State Problems 

The discussion of global functions In the two earlier sections 
pertained to elastostatlc analysis of the localized stress zones. Here 
some remarks on steady-state dynamic effects are made. The one-dimensional 
finite-element method for generating global function data bases can be 
modified for steady-state Inertial effects by Including kinetic energy In 
the problem formulation. Instead of Eqs. (4) and (5)» those equations 
become* respectively: 

CK 1 3 {Q> - yCK 2 3{Q} + T 2 [K 3 3{QJ + w 2 [M]{Q) = 0 (6) 

[Kj(m)]{QJ - Y[K 2 (m)]{Q) + Y 2 [K 3 (m)]{Q} + « 2 [M]{Q) = 0 (7) 

where w Is the steady-state forcing frequency. The derivations of these 
equations are given In Refs. [11*143. 

With these global functions* It Is possible to study elastic wave 
scattering In prismatic structures by discontinuities during vibration 
or by some other steady-state dynamic Input. GLFEM analysis of this 
type of prismatic structures will be similar to problems of elastic wave 
scattering by an object embedded In an Infinite medlun or soil -structure 
1 nteractl on. 


Effects of Initial Stress 

Using the same methodology* prismatic structures under Initial 
stress may also be analyzed. In this case* the global functions must 
Include the prestressing effect. One-dimensional finite-element 
analysis of wave propagation In laminated composite plates and cylinders 
under Initial stress have been explored* see Refs. [15*163. It Is a 
straightforward task to adapt these formulations to generate an 
elgenproblem for the global functions for a prismatic structure under 
Initial stress. Also, no conceptual difficulties are seen In an 
extension to three-dimensional prismatic structures under prestress. 


Concl udl ng Remarks 

Considerable discussion has been devoted to the strategies of GLFEM 
analyses of prismatic structures with localized stress regions and other 
discontinuities. The role of the global functions has been clearly 
outlined and their mathematical requirements Indicated. The method 
for deriving these global functions for prismatic structures, whose 
cross-sectional geometries are complicated by laminated construction* 
has been discussed. From the discussion of GLFEM analysis strategy. 

It should be clear that GLFEM Is feasible and effective. Considerable 
economy of computational efforts over a strictly FEM approach should 
be real Ized. 
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Basic global-local mesh configurations. 
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. Elastic wave scattering by axi symmetric inclusion embedded in 
an infinite isotropic medium. 
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TRACTION-FREE 



Figure 3. Soil -structure interaction problem. 
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SCARF JOINT OF ALUMINUM AND COMPOSITE CYLINDERS 


Figure 4. Two-dimensional prismatic structures. 
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in a composite plate. 
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EXAMPLE OF A 

THREE-DIMENSIONAL PRISMATIC STRUCTURE 



Figure 6. Three-dimensional prismatic structures. 
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